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We investigated the role of oxygen vacancy in n-type interface of LaAlOs (LAO) overlayer on SrTiOs (STO) 
(001) by carrying out density-functional-theory calculations. Comparing the total energies of the configurations 
with one vacancy in varying locations we found that oxygen vacancies favor to appear first in LAO surface. 
These oxygen vacancies in the surface generate a two-dimensional distribution of carriers at the interface, re- 
sulting in band bending at the interface in STO side. Dependent on the concentration of oxygen vacancies in 
LAO surface, the induced carrier charge at the interface partially or completely compensates the polar electric 
field in LAO. Moreover, the electronic properties of oxygen vacancies in STO are also presented. Every oxygen 
vacancy in STO generates two electron carriers, but this carrier charge has no effect on screening polar field in 
LAO. Band structures at the interface dependent on the concentrations of oxygen vacancies are presented and 
compared with experimental results. 

PACS numbers: 73.20.-r, 79.60. Jv, 77.22.Ej, 73.21. -b 



I. INTRODUCTION 

Since Ohtomo and Hwang reported the existence of a high 
mobility electron gas at the n-type (LaO)+/(Ti02)^ interface 
between two band-gap insulators LAO and STO\ many novel 
properties related to the electron gas at the interface, such as 
insulator-metal transition, superconductivity, and ferromag- 
netism were found^i^i^i^. The mechanism of conductivity 
and the dimensionality of the induced carrier at the interface 
were intensively studied. In the n-type interface structure a 
polar electrical field along LAO [001] direction is arisen from 
the alternating stack of (LaO)+ and (A102)~ charged lay- 
ers. Charge reconstruction at the interface was proposed as a 
way to avoid diverging electrostatic potential with thick LAO, 
which compensates the polar electric field in LAO and induces 
electron gas at the interfacai*^. The densities and distribu- 
tions of the induced carriers charge were measured in various 
experiment a^iVi^i^Qi^^^^i^^i^'^1^^1^^^'^ . Two different origins of 
the induced carriers charge, i.e. intrinsic and extrinsic doping, 
have been proposed. Intrinsic charge doping at the interface 
mainly happens in stoichiometrical structure without oxygen 
vacancy^iLiSii^. With thick LAO overlayer the charge trans- 
fers from valence band of LAO surface to conduction band 
of STO at the interface, leading to a metallic interfacdi2i22i^. 
While extrinsic charge doping was found to be ascribed to 
oxygen vacancies in many experiments 

Noticeably, the conducting properties at intrinsic and ex- 
trinsic doped interfaces measured in the experiments are 
very different. In the ideal interface of thin LAO overlayer 
without oxygen vacancy the carrier densities measured by 
Hall effect and x-ray photoelectron spectroscopy (XPS) are 
even less than 0.1 electron per two-dimensional unit cell (2- 
d.u.c.)2ii2iii. Recent density-functional-theory (DFT) calcu- 
lations showed that ionic polar distortions partially compen- 



sate polar electric field in LAO, reducing carrier density at 
the interfacei2i21i22i22i2i2^. For thin LAO overlayer the carrier 
charge was found in XPS experiment to distribute in a few 
layers of STO at the interface—. Latest DFT study found that 
for less than 7 layers of LAO the carrier charge distributes in 
3 layers of STO but for thicker LAO the carrier could accu- 
mulate in deep STO layers about 3 nm away from the inter- 
face. However, with oxygen vacancies in the interface struc- 
ture the carrier densities measured in experiments vary from 
0.5 to three orders of magnitude e/2-d.u.c. dependent on the 
oxygen pressure in preparing proces&ii^ii^iil. The experiment 
in Ref.l 12] explicitly verified that high carrier density is as- 
cribed to high concentration of oxygen vacancies in STO sub- 
strate which were generated under lower oxygen pressure dur- 
ing preparing the sample. 

Besides the carrier density, another issue which people 
have been debating extensively is the distribution of carrier 
charge induced by oxygen vacancy. In some sample pro- 
duced in high oxygen pressure two-dimensional distributions 
were observed^'^'^^'^^, but in others produced in low oxy- 
gen pressure obvious 3 -dimensional distribution in STO were 
observediSii^. These experimental hints indicated the carrier 
distribution is related to the concentration of oxygen vacan- 
cies in some way. Up to now, people have realized that oxy- 
gen vacancy play an important role in n-type LAO/STO in- 
terface, however the locations of oxygen vacancies and the 
dependence of two-dimensional and three-dimensional carrier 
distributions on the locations of oxygen vacancies are still not 
clear. 

In this paper we studied the n-type interface with oxygen 
vacancies by using density-functional-theory (DFT) calcula- 
tions. It should be noted that the interface composed of sand- 
wich structure was also investigated in experiment but the 
structure of LAO overlayer on STO substrate has attracted the 
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most of attentions because of its abundant conducting proper- 
ties. In this paper we focused on the structure of LAO over- 
layer on STO. We calculated and compared the total energies 
of the configurations with varying locations and concentra- 
tions of oxygen vacancies. We found that the locations of 
oxygen vacancies are related to the concentration. From the 
calculations of electronic structures we also unveiled that the 
distribution of carriers is strongly dependent on the locations 
of oxygen vacancies. Based on the results from our calcu- 
lations, energy band structures at the interface and localized 
surface conductivity were discussed. 
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FIG. 1: (Color online) Atomic structures of (2x2) supercells (a) 
without oxygen vacancy, and (b) with oxygen vacancy in LAO sur- 
face. Oxygen vacancy is denoted by black ball. 



II. COMPUTATIONAL METHOD 

We carried out DFT calculations by using the Vienna ah 
initio Simulation Package (VASP)^^ within a generalized gra- 
dient approximation^^ together with the projector augmented 
wave pseudopotentials^^'^^ and the cut-off energy of 400 eV 
for the plane wave basis. We modeled the LAO/STO interface 
by a slab consisting of 2 to 7 LAO layers on top of STO(OOl) 
substrate and a vacuum region of 14 A along the c-axis in a 
supercell geometry. Dipole corrections were used to correct 
the errors of electrostatic potential, forces, and total energy 
caused by periodic boundary condition^. To simulate vari- 
ous concentrations of oxygen vacancies the in-plane unit-cells 
were taken as (2x2), (3x2), and (3x3). F-centered (5x5), 
(3x5), and (3x3) k-point meshes were used to sample the 
Bullion zone, respectively. Fig.[T](a) and (b) present relaxed 
(2x2) supercells without oxygen vacancy and with oxygen 
vacancy in LAO surface, respectively. The in-plane lattice 
constant of the slab was constrained at the calculated equi- 
librium lattice constant a = 3.942 A of the STO substrate. 



All coordinates of atomic positions were fully relaxed with 
forces less than 0.02eV/A except for the atoms in the bottom 
two layers of STO, which were fixed in their bulk positions. 



III. LOCATION OF OXYGEN VACANCIES OF LOW 
CONCENTRATION 

In this paper we divide the concentrations of oxygen va- 
cancies in the interface structure into two regions, i.e. low 
concentration and high concentration. In the former the con- 
centrations of oxygen vacancies is less than or equal to 1/4 
per 2-d.u.c. and the density of induced carriers is less than or 
equal to 0.5 electron per 2-d.u.c.. While in the latter the con- 
centrations of oxygen vacancies is larger than 1/4 vacancy per 
2-d.u.c. and the density of induced carriers is larger than 0.5 
electron per 2-d.u.c.. To simulate low concentration of oxygen 
vacancy in the interface structure only one vacancy is involved 
in the supercells. While for higher concentration more vacan- 
cies are involved in the supercells. For convenient we use ny 
and nc to denote the concentration of oxygen vacancies and 
carrier density in 2-d.u.c., respectively. 



A. Binding energies of oxygen atom 

Oxygen vacancies were explicitly confirmed to exist in 
STO substrate of the interface structure produced under low 
oxygen pressur e However, up to now there has been no 
report about whether oxygen vacancy exits in LAO. To inves- 
tigate the location of oxygen vacancy we calculated the total 
energies of oxygen vacancy in varying locations under vari- 
ous concentrations. We used (3 x 2) and (2 x 2) supercells with 
one vacancy as the representatives of ny less than and equal to 
1/4 per 2-d.u.c., respectively. To compare them in same plot 
we calculated the binding energy of one oxygen atom corre- 
sponding to the vacancy defined as: 



^5 = (^y + l/2^oj-^o, 

in which Ey, Eq^, and £^0 are total energies of the system 
with oxygen vacancy, oxygen molecule and ideal system with- 
out vacancy, respectively. Figure [2] presents the binding ener- 
gies of one oxygen atom corresponding to the vacancy at vary- 
ing locations in (3x2) and (2x2) supercell. Oxygen atom in 
LAO surface has the lowest binding energy, which is far less 
than that in STO substrate. This indicates that oxygen vacancy 
is most easily formed in LAO surface rather than in STO sub- 
strate. Moreover the binding energy of oxygen atom in LAO 
surface decreases with the thickness of LAO. Such results are 
strongly related to polar electric field in LAO. In the ideal in- 
terface structure the polar field shifts the energy bands of LAO 
toward higher energy layer by layer, substantially raising the 
total energy. While oxygen vacancy in LAO dopes electron 
carriers at the interface, which partially or completely screen 
the polar field in LAO, remarkably lowering the total energy. 
In contrast to the vacancy in the surface, the vacancy inside 
LAO just screens part of polar field from the layer with the 
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vacancy to the interface. Therefore corresponding configura- 
tion has higher total energy, as illustrated in Fig. [21 
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Location of oxygen vacancy 

FIG. 2: (Color online) Binding energies of oxygen atom correspond- 
ing to one vacancy at varying locations in (3 x 2) supercell consisting 
of STO substrate and 4 LAO layers and (2 x 2) supercells consisting 
of STO substrate and 2 or 3 LAO layers. 



B. Formation energy of oxygen vacancy 
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FIG. 3: (Color online) Formation energies of one oxygen vacancy in 
LAO surface as a function of oxygen chemical potential. The thick- 
nesses of LAO vary from 2 to 5 layers for (2x2) supercells, and 
only 4 layers of LAO for (3x2) supercell. Formation energy of ideal 
structure without vacancy is taken as zero. 

To investigate dependence of the formation of oxygen va- 
cancy in the surface on the thickness of LAO, we calculated 
the formation energies of one vacancy in the surface of (2 x 2) 
supercells consisting of 2 to 5 layers of LAO according to the 
formula 

Ej = Ey — {Eo — /io), 

in which /io is chemical potential of oxygen atom, Ey, and 
£^0 are total energies of the system with one oxygen vacancy 
in the surface and ideal system without vacancy, respectively. 



As shown in Fig.[3]the formation energies of (2x2) supercells 
with one vacancy in LAO surface decrease with the thick- 
ness of LAO overlayer and trend to converge for those of 
thickness of LAO more than 4 layers. As aforementioned, 
in the ideal structure the electrostatic energy induced by po- 
lar field increases with the thickness of LAO. Following elec- 
tronic calculations show that for (2x2) supercells with one 
oxygen vacancy in the surface the density of doping charge at 
the interface is 0.5 e/2-d.u.c., which exactly compensates po- 
lar electric field. This implies that oxygen vacancy in the sur- 
face of thicker LAO overlayer eliminates more electrostatic 
energy than that of thinner LAO overlayer. While once the 
LAO overlayer exceeds 4 layers, as demonstrated in previ- 
ous DFT studies, in ideal structure without vacancy the intrin- 
sic doping charge at the interface partially screens the polar 
field, preventing the increasing of electrostatic energy with the 
thickness of LAO^W^ Therefor, the enerRv difference be- 
tween vacant and ideal structures does not increase with thick- 
ness of LAO. So we can predict that under same preparing 
condition oxygen vacancy can be formed more easily in the 
sample with thicker LAO film, while this tendency does not 
strengthen any more once LAO exceeds 4 layers. Moreover, 
as shown in Fig.[3l for two unitcells (2x2) and (3x2) which 
have the same 4 layers of LAO the lower concentration of va- 
cancy is more stable and therefor more easily be formed. As 
well known, oxygen pressure is a dominant factor of the for- 
mation of oxygen vacancy in LAO/STO interface structure. 
In Fig. [3] at oxygen rich limit the ideal structure has the lowest 
formation energy. This is consistent with the experimental re- 
sults that vacancy-free structures were fabricated under high 
oxygen pressure^iiSii. While in the middle or near the poor 
limit vacant structures have lower formation energies. 



IV. ELECTRONIC STRUCTURES OF OXYGEN 
VACANCIES OF LOW CONCENTRATION 



A. Oxygen- vacancy state in LAO surface and induced carrier 

Above binding-energy calculations have shown that in the 
configurations of ny <=l/4 per 2-d.u.c. oxygen vacancy lies 
in LAO surface. To investigate the electronic properties of 
oxygen vacancy in the surface we calculated (2x2), (3x2), 
(3x3), and (4x4) supercells with one vacancy, respectively. 
As a representative of ny <l/4 per 2-d.u.c., the electronic 
structure and local atomic structure of (3x3) supercell are 
presented in Fig. (H One can see in Fig. [4] (a) and (b) that 
the empty state of oxygen vacancy lies a little bit below the 
conduction band minimum (CBM) of LAO and is localized in 
the surface layer. Obviously, the level of oxygen vacancy is 
higher than CBM of STO. Charge transfer from the state of 
oxygen vacancy in the surface to conduction band of STO at 
the interface lowers the total energy of the system. Our cal- 
culation showed that 2 electrons accumulate at the interface, 
i.e. nc= 0.22 e/2-d.u.c., which is the twice of ny. The same 
relation between ric and ny were also found in other configu- 
rations of ny <=l/4 per 2-d.u.c.. 
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FIG. 4: (Color online) (a) Layer-projected DOS for the (3x3) super- 
cell with one oxygen vacancy in LAO surface, (b) Spacial distribu- 
tion of oxygen- vacancy state in the (3 x 3) surface, (c) Side view of 
polar distortions of cations and anions in upper 3 layers of LAO. 



B. Lattice polarization dependence on the concentration of 
vacancy in LAO surface 

For the configurations of ny <l/4 per 2-d.u.c. the carrier 
density at the interface is less than 0.5 e/2-d.u.c.. The polar 
electrical field in LAO overlayer is screened partly and then 
the potential gradient in LAO overlayer still remains. In Fig.lH 
(a) one can see that the valence and conduction bands shift 
towards higher energy layer-by-layer from the interface to the 
surface. The polar distortions in LAO induced by residual 
polar electrical field are shown in Fig.|4](c). Summarizing for 
configurations of ny <=l/4 we found that with increasing 
the concentration of oxygen vacancy from (4 x 4) to (2 x 2) the 
density of carriers at the interface increases till 0.5e/2-d.u.c 
and the potential gradient and the polar distortions reduce till 
vanish. 

For the configurations of ny <l/4, the thickness of LAO is 
another factor which affects the density of carriers at the in- 
terface. For a given concentration of oxygen vacancies, with 
increasing the thickness of LAO the VBM of LAO shifts up- 
wards until it touches CBM of STO. This critical thickness 
of LAO, denoted as tc, is dependent on the concentration of 
vacancies in the surface ny. Once the thickness of LAO over- 
layer overruns tc the VBM of LAO does not shift up anymore, 
instead, more charge transfers from valence band of LAO sur- 
face to the interface. According to the electrostatics we give 
an expression of carrier density nc as a function of ny and t 
(thickness of LAO overlayer) 



where A=1.97, to=0.053 obtained from our previous DFT 
study. The scale of t and to is unit-cell layer. For the con- 
figuration of (3x3) with one vacancy in the surface the criti- 
cal thickness equals 7 u.c in our DFT calculation and 7.1 u.c. 
from above formula. 

Figure[5](a) presents the layer-projected DOS s of (2x2) su- 
percell with one vacancy in LAO surface. The state of vacancy 
in LAO surface Hes a little bit below the CBM of LAO and is 
localized in the surface layer. Calculated density of carriers at 
the interface is equal to 0.5 e/2-d.u.c., which is exactly enough 
to completely screen the polar electric field in LAO. One can 
see this point from straight- aligned VBMs and CBMs of LAO 
layers in Fig. [5] (a). Similar to what is shown in Fig. \T\ (b), 
oxygen vacancy in the surface induces strong tilt distortion 
which even extend into a few layers of STO, but average polar 
displacement in each LAO layer is zero. It means no net elec- 
trostatic field in LAO. For comparison Fig.[5](b) presents the 
layer-projected DOSs of (2x2) supercell with one vacancy in 
STO substrate. Oxygen vacancy in STO produces a shallow 
energy level in the CBM of STO and generates two electron 
carriers. However, strong potential gradient appears in LAO 
overlayer because no charge transfers from the LAO surface 
to the interface. We can conclude that carriers generated by 
oxygen vacancies in STO hardly have any effect on screening 
the polar electric field in LAO overlayer. 
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FIG. 5: Layer-projected DOS for the (2x2) supercell (a) with one 
oxygen vacancy in LAO surface, and (b) with one oxygen vacancy in 
STO-1 layer. 



DISTRIBUTION OF OXYGEN VACANCIES OF HIGH 
CONCENTRATION 



A. Binding energies 
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We modeled the configurations of ny >l/4 per 2-d.u.c. by 
(2x2) supercell with two oxygen vacancies. As we pointed 
out above that when ny <=l/4 per 2-d.u.c. oxygen vacancy 
favors to lie in LAO surface. Considering this point in the cal- 
culations of (2 X 2) supercell with two oxygen vacancies one is 
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FIG. 6: (Color online) Binding energies of one oxygen atom at vary- 
ing locations in (2x2) supercell consisting of 4 u.c. layers of LAO 
under the condition of one oxygen vacancy fixed in LAO surface. 
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always fixed in LAO surface and the others varies from LAO 
surface to STO subxtrate. Fig. [6] illustrates binding energies 
of the oxygen atom corresponds to the second vacancy. Al- 
though oxygen atom in surface still has the lowest binding 
energy, in contrast to the cases of ny <=l/4 per 2-d.u.c. the 
energy differences with other configurations decrease largely. 
The energy difference between the configurations of the sec- 
ond vacancy in LAO surface and in STO substrate is about 
0.2 eV, which is far less than that for the first vacancy. This 
implies that second vacancy has large probability to appear 
in STO. Following electronic structure calculations show that 
for the configurations of ny >l/4 per 2-d.u.c. the second va- 
cancy does not contribute any energy gain by screening the 
electrostatic field in LAO because the first vacancy already 
completely screens it. 
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FIG. 7: Layer-projected DOSs for (2x2) supercells (a) with two oxy- 
gen vacancies in LAO surface, and (b) with one vacancy in LAO 
surface and one in LAO-3 layer. 



FIG. 8: Layer-projected DOSs for the (2x2) supercells (a) with one 
oxygen vacancy in LAO surface and one in STO, and (b) with one 
oxygen vacancy in LAO surface and two in STO. 



B. Electronic structure of oxygen vacancies of high 
concentration 

Figure|7](a) presents the layer-projected DOS of (2x2) su- 
percell with two oxygen vacancies in surface AIO2 layer. The 
state of vacancies is widened due to the interaction of two 
vacancies, and only half is occupied. Two electrons in this 
state transfer to the interface. Another typical configuration is 
shown in Fig.jTjb), in which one vacancy lies in surface AIO2 
layer and the other in the (LAO-3)-A102 layer. The state of 
oxygen vacancy in LAO-3 layer extends into the surface layer 
and is fully occupied, while the state of oxygen vacancy in the 
surface layer is empty. Our calculated amount of carriers at 
the interface is still 2. In these two configurations the former 
has a metallic surface, while the latter has an insulating sur- 
face. Both dope two electron carrier at the interface, although 
two oxygen vacancies exist in LAO. Several other configura- 
tions with two oxygen vacancies in LAO also show the same 
carrier density 0.5e/2-d.u.c at the interface. This indicates that 
oxygen vacancies in LAO can contribute at most 0.5 electron 
carrier per 2-d.u.c. at the interface, and extra vacancies apart 
from 1/4 per 2-d.u.c. have no contribution to the conductivity 
at the interface. It can be understood from the sight of polar 
electric field in LAO. Charge transfer of 0.5 e/2-d.u.c from the 
LAO surface to the interface completely compensate the polar 
field in LAO, lowering the total energy. While Charge trans- 
fer of more than 0.5 e/2-d.u.c would set up an inverse electric 
field in LAO, which could raise the total energy. 

Since a great amount of oxygen vacancies were found in 
STO substrate in the samples produced under lower oxygen 
pressure, the contribution of vacancies in STO to density and 
distribution of the carriers should be investigated. Fig. Oa) 
shows the Layer-projected DOSs for the (2x2) supercell with 
one oxygen vacancy in LAO surface and one in STO. The car- 
rier density is 1.0 e/2-d.u.c. in our calculation. With one more 
vacancy in STO, as shown in Fig. Ob), the carrier density in- 
creases 0.5 e/2-d.u.c. more. In contrast with the upper limit 
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contribution of the vacancies in LAO every vacancy in STO 
generates two electron carriers in STO. Considering the va- 
cancies distributing uniformly in whole STO substrate under 
low oxygen pressure in experiments, a 3-dimensional distri- 
bution of carriers can be formed. In respect of band structure 
no potential gradient appears in LAO, while valence and con- 
duction bands of STO have evident bending at the interface. 



VI. DISCUSSION 

A. Band diagrams 

According to above calculated electronic structures we plot 
schematic band diagrams of the interface. Since oxygen va- 
cancies in STO do not cause the band bending at the interface, 
only the band structure with various concentrations of oxy- 
gen vacancies in LAO are plotted. Fig. [9] presents the band 
structures of oxygen vacancies in or near LAO surface. At 
the interface in STO side strong band bending happens in all 
configurations. The reason is that net positive-charged LAO 
overlayer, produced by electrons transfer from LAO surface 
to STO, generates an attractive potential to electrons in STO, 
lowing the potential at the interface in STO side. In Fig. [9] 
(a) ny <l/4 per 2-d.u.c. and nc < 0.5 e/2-d.u.c.. The polar 
electrical in LAO is partially screened and the residual field 
results in bands sloping in LAO. In Fig. [9](b) ny =1/4 per 
2-d.u.c. and nc = 0.5 e/2-d.u.c.. The polar field is com- 
pletely screened, leading to flat bands in LAO. In contrast with 
ny <l/4 per 2-d.u.c. in the case of ny =1/4 per 2-d.u.c. more 
charge transfers from the surface to the interface, leading to 
larger band bending at the interface in STO side. In Fig.[9](c) 
ny >l/4 per 2-d.u.c. but still nc = 0.5 e/2-d.u.c.. The band 
structure is the same as the case of ny =1/4 but the surface is 
metallic. In Fig. [9] (d) vacancies are in surface and the layer 
below surface and the surface is insulating. In previous XPS 
experiment a flat valence band was observedi^. According to 
above analysis one can easily realize that such a band structure 
is arisen from the oxygen vacancies of ny >=l/4 per 2-d.u.c. 
in LAO surface. 



B. Charge carriers distribution: 2-dimensional or 
3-dimensional 



To clarify the distribution of the carriers induced by oxygen 
vacancy we carried out the calculations of electronic structure 
for (LAO)3/(STO)i7-(2xl) supercefl by employing 17 STO 
layer. With one oxygen vacancy in LAO surface the integrated 
carrier density in STO is 0.5 e/2-d.u.c.. While the carriers 
in STO consist of two separate components, as illustrated in 
Fig.[T3a). The interface component contains the main part of 
carriers, and the extended part with a broad peak locates near 
the 14th layer, i.e. about 5 nm, from the interface. The po- 
tential profile in the STO side illustrated in Fig. [9] resembles 
the case of the inversion layer in metal-oxide- semiconductor 
field-effect transistor and semiconductor heterostructure^. 
Within this triangle-like potential, the carriers in the lowest 
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FIG. 9: Band diagram of the interface of LAO overlayer on STO with 
oxygen vacancies in LAO surface for (a) ny < 1/4, (b) ny = 1/4, 
(c) ny > 1/4, (d) ny > 1/4 but vacancies in both surface layer and 
the layer below surface. 
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FIG. 10: (a) Layer-resolved carrier density in STO for 
(LAO)3/(STO)i7-(2x 1) with one oxygen vacancy in LAO surface, 
(b) Charge density plot of the induced carrier. 



bound state which can be approximately described by the Airy 
function accumulate at a few nm away from the interface^. 
However, the interface state is originated from different mech- 
anism. To understand the nature of two different components 
of the carriers, we plotted the charge density of the induced 
carrier in the STO side in Fig. [TOlb). The interface compo- 
nent consists mostly of the Ti dxy orbitals, but the extended 
component has contributions from all the Ti t2g orbitals, i.e., 
dxy, dyz, and dzx states. The character at the interface stems 
from a strong compressive distortion of the TiOe octahedron 
at the interface, which produces a strong tetragonal field, low- 
ering of the dxy state. Similar character of carriers was found 
in the ideal interface of LAO/STO without vacancy, and the 
mechanism was detailedly analyzed in our previous work^. 

From Fig. [TOlone can see that most of carriers accumulate 
at the interface in four layers of STO although there is ex- 
tended state. Such a two-dimensional distribution of carri- 
ers induced by oxygen vacancy in LAO surface is consistent 
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with the two-dimensional electron gas measured in LAO/STO 
samples produced under higher oxygen pressure While 
under low oxygen pressure oxygen vacancies can be formed 
uniformly in STO besides LAO surface. As well known, oxy- 
gen vacancies in STO generate a shallow level in conduction 
band and Fermi level lies in the CBM of STO. Considering 
band-bending at the interface the carriers first fill in the po- 
tential valley at the interface and then in STO inside. This 
implies that large amount of oxygen vacancies in LAO/STO 
system would generate a three dimensional distribution of car- 
riers in STO but with higher carrier density near the inter- 
face. Above picture about carrier distribution was measured 
in previous experiment by using conducting atomic forcer mi- 
croscopy (AFM), in which the carrier density was found to 
decays exponential near the interface and extends a few fim in 
STOi^. 



C. Localized surface conductivity 

Another finding in our calculations which we like to em- 
phasize is the surface conductivity in the system of LAO/STO. 
As shown in Fig.[7ta), oxygen vacancies of more than 1/4 per 
2-d.u.c. generated a metallic surface. In contrast to nearly uni- 
form distribution of the carriers at the interface, the carriers in 
the surface is relatively localized. In the case of Fig. Ha), 
the states of two close vacancies overlap. A part of electrons 
in the states transfer to the interface and the left form local- 
ized surface carriers. This localized surface carriers may be 
responsible for the few-nanometer- size conducting islands in 
the surface of LAO/STO found in recent experiment^. 



VII. SUMMARY 

We investigated the n-type LAO overlay er on STO(OOl) 
substrate with oxygen vacancies by using DFT calculation. 
We found that oxygen vacancies favor to appear first in LAO 
surface, which generate a two-dimensional carriers at the in- 
terface. The density of carriers induced by vacancies in LAO 
surface has an upper limit 0.5 e/2-d.u.c.. For the concentra- 
tion of vacancies in LAO surface less than 1/4 per 2-d.u.c., 
the density of induced carriers at the interface is less than 0.5 
e/2-d.u.c. and the energy bands in LAO slope up from inter- 
face to surface. While for that of equal to or more than 1/4 
per 2-d.u.c. the density of induced carriers is equal to 0.5 e/2- 
d.u.c. and no slope of the energy bands occurs in LAO. We 
found that for the case of the concentration of oxygen vacan- 
cies more than 1/4 per 2-d.u.c. the surface presents metallicity. 
We also investigated the role of oxygen vacancy in STO. We 
found that every oxygen vacancy in STO generates two elec- 
tron carrier, but this carrier charge has no effect on screen- 
ing the polar electric field in LAO. We predict that when a 
large amount of oxygen vacancies present in LAO/STO sys- 
tem the carriers in STO show a three-dimensional distribution 
with higher density at the interface. 
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